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I. INTRODUCTION 

Superconducting quantum bits (qubits) [T] are a lead- 
ing candidate for a solid-state quantum computer. How- 
ever, while coherence times are continually increasing, it 
remains necessary to study how to maximize coherence 
while accessing the large Hilbert space required by key 
applications in quantum information processing. Exam- 
ples include the rapid controUed-phase gate using auxil- 
iary states [IHS] and the general framework to improve 
quantum logic gate synthesis using multilevel systems 
[6j. An emerging pattern is that resources outside of 
the traditional qubit states can lead to improved control 
sequences with reductions in total time or complexity. 

Superconducting phase and transmon circuits are a 
natural candidate to explore operations outside of the 
qubit subspace, as these systems are in fact weakly anhar- 
monic oscillators with many levels. Control of multiple 
levels in these devices has been demonstrated experimen- 
tally [THTU] and explored theoretically [TTHIS] • Notably, a 
theoretical method [16] incorporating multiple levels has 
led to improvements in qubit logic operations [T71 [TB] . 

Superconducting resonators can also be controlled at 
the Fock state level. By coupling such resonators to 
an auxiliary nonlinear system, recent experiments have 
created Fock states [TU], observed their decay [5D], and 
demonstrated the synthesis of arbitrary superpositions of 
Fock states [H]. This last experiment used the protocol 
of Law and Eberly |22j with linear coupling of a phase 
qubit to a coplanar waveguide resonator. A recent theo- 
retical work 2^1 extended this approach to the synthesis 
of entangled states of two (and possibly more) resonators. 
Subsequently, an experimental synthesis [24] of a "high" 
NOON state [5S] was accomplished using an alternative 
procedure 

While great progress has thus been made in the con- 
trol of superconducting resonators, these works leave 
open the question of whether the larger state space of 
the resonator can be used to process quantum informa- 
tion. While there are certainly caveats to the question of 
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"qubit or oscillator?" (see, e.g. [H]), there is an estab- 
lished body of work demonstrating that quantum sys- 
tems with multiple states, known as qudits (for d-level 
systems), can be as useful as qubits. Using the lowest 
d levels of a harmonic oscillator would thus be a poten- 
tial alternative to qubits. For superconducting circuits 
in particular, it is clear that resonators can be fabricated 
with much greater precision and coherence, and thus a 
central question is how to compute using the additional 
resources present in harmonic oscillator modes. 

An important step in that direction was taken by Ja- 
cobs [2H] , who showed that linear coupling of an oscillator 
to an auxiliary qubit was sufficient to approximate any 
desired evolution of the oscillator. This was based on 
the general Lie algebraic result by Lloyd et al. [29], but 
applying this to quantum logic on a discrete set of Fock 
states would require significant overhead in complexity 
(to synthesize the desired interactions). Inducing a non- 
linearity perturbatively | 30j is another route to unitary 
control of the oscillator, although this may require some 
compromise in timescales (to stay within the perturba- 
tive limits). A scheme of this sort, appropriate to atomic 
cavity-QED or ion trap systems, was proposed by Santos 
|31| and serves as a primary inspiration for our proposal. 
Here we present a detailed analysis of a circuit-QED ap- 
proach in which a three-level system, such as a phase or 
transmon qubit, is used as an auxiliary to enable arbi- 
trary unitary control of a superconducting resonator. 

In this work we combine two experimentally demon- 
strated interactions to propose a simple procedure to 
perform an arbitrary rotation between Fock states, and 
by composition an arbitrary unitary operation on the 
Fock states. The basic idea is shown in Fig. [l] We 
devise a control sequence to selectively move two states 
of the oscillator to auxiliary levels. Here a rotation or 
swap S is performed between these two auxiliary lev- 
els, and finally these levels are returned to the original 
oscillator states. For convenience, we will call the first 
step an encoding operation C/cncodo and the final step a 
decoding operation J/dccodo, so that the net rotation is 

C^l,2 = CAiecode'S'i/encodo- 

The first ingredient in our proposal is a quasi- 
dispersive interaction between a qubit and the resonator, 
to allow for number-state-dependent rotations of the 
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FIG. 1: General approach to quantum logic using harmonic 
oscillator states. Encoding and decoding operations transfer 
two oscillator states (here n = 1 and n = 2) to a pair of 
auxiliary states which can be swapped by S. Here the total 
single-qudit rotation is t/1,2 = C^decode'S'C^ncode 



qubit. This was first seen spectroscopically by Schus- 
ter et al. |32j . and more recently used to perform a 
non-demolition measurement of a resonator memory by 
Johnson et al. [33j . We propose to use this interaction 
to selectively address the Fock states of interest as part 
of the encoding and decoding operations. This approach 
was previously used in the entangled state synthesis al- 
gorithm |23j . 

The second ingredient is a resonant swapping inter- 
action between the resonator and the higher levels of a 
superconducting phase or transmon qubit. This was used 
in the NOON-state synthesis experiment and effects 
the qudit rotation by swapping the auxiliary levels, as 
shown in Fig. 1. 

This paper is organized as follows. In Section II we 
briefly review existing qudit theory and outline how our 
scheme can be used for qudit logic operations. In Sec- 
tion III, the basic system of a three-level system coupled 
to an oscillator is presented and analyzed. In Section 
IV the time-dependent control sequences are described 
and verified by numerical simulations. In Section V we 
show how this can be extended to a two-qudit logic gate. 
In Section VI we analyze the effects of decoherence and 
discuss resonator measurement. Finally, we conclude in 
Section VII with a discussion of open topics for study. 



II. QUDIT LOGIC 

Multilevel quantum logic has been explored as an al- 
ternative to the traditional qubit constructions by many 
authors [34H37] . We follow the discussion by Brennen et 
al. |38) . They show, using the QR decomposition from 
linear algebra, that arbitrary single-qudit unitaries can 
be constructed from a family of two-component rotations 



where we have defined two operations in the qudit sub- 
space |fc)}: 



exp 



and 



7^|.fc(0) = cxp 



-\k){k\) 



(2) 



(3) 



In addition to these single-qudit rotations, we also need 
a two-qudit operation, to generalize the controUed-NOT 
gate commonly used in qubit circuits. We shall synthe- 
size the controUed-phase gate 



U,A0)^exp{~t9\j,k){j,k\). 



(4) 



where \ j, k) = |j) ®\k) is the state in which the first qudit 
is in state \i) and the second qudit is in state \k) . This 
gate set is sufficient to perform an arbitrary two-qudit 
unitary operation |38| . and by extension to multiple qu- 
dits, universal quantum computation [35]. Explicit con- 
structions for circuit synthesis can be found in [33 HO] ■ 

In the implementation we will present shortly, the ro- 
tations will be between neighboring oscillator states j 
and k = i + 1. That is, we will construct C/jj+i from 
the sequence ?7dccodc'5'C4ncodc, as illustrated in Fig. 1, 
where S performs a swapping interaction between the 
amplitudes for states |j) and \ j + 1). Note that this limi- 
tation to neighboring oscillator states does not present a 
true obstacle to general qudit logic. As shown in Lemma 
II. 1 of j38j . the important requirement is that there is a 
connected coupling graph between the qudit states. Ro- 
tations between neighboring states leads to a linear cou- 
pling graph: 



0olo2o---od-l. 



(5) 



(1) 



Finally, the single-qudit phase rotations can be 
performed as in current experiments [21 13] , by short de- 
tuning pulses that can be incorporated in the single-qudit 
rotations. 



III. SUPERCONDUCTING IMPLEMENTATION 

We extend the framework of [53], in which two super- 
conducting resonators are coupled by a tunable circuit, 
as shown in Fig. [2j Letting a and be the ladder oper- 
ators for resonator A, h and for resonator i?, we model 
the system as 

T-L = Ha + hujaa^a + hga {aa^ + aV^) 

+ Hb + fiLUbb'^b + hgb {ba^ + b^^a^) 

+ hgab {afa^ + c7-^af) ■ (6) 

Here Ha and Hb are the single-qubit Hamiltonians for 
the auxiliary, and a±. and are the corresponding rais- 
ing and lowering operators (see below). We will as- 
sume the coupling gab between the two auxiliaries can 
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FIG. 2: Schematic superconducting circuit to implement sin- 
gle and two-qudit operations on resonators A and B. Each 
resonator is coupled to an auxiliary qutrit {A or B), which 
are themselves coupled to each other. Control lines allow ma- 
nipulation of the qutrits. 

be turned on and off at will, using the tunable coupling 
circuits recently demonstrated [HI , and that the aux- 
iliaries can be controlled by microwave and flux pulses. 



(a) - (b) 
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FIG. 3: (a) Energy level diagram for a three-level artificial 
atom coupled to a resonator. The artificial atom has level 
spacings loqi and a;i2 < cjqi , while the resonator has frequency 
ujr (b) Approximate energy level diagram in the dispersive 
regime, in which the second excited state of the atom has been 
eliminated. The dressed states have number-state-dependent 
level spacings ojq"' . 



A. Single Resonator Model 

We begin by focusing on a single qubit-resonator sys- 
tem (i.e. Aor B), described by the following Hamiltonian 

y. = Hq + hLUro'a + hg (0(7+ + aV_) , (7) 

where the auxiliary quantum system is taken as a three- 
level qutrit. 



and 




nQ=h\ cjoi 
UJ02 



1 
a = \ A 




(8) 



(9) 



with (T+ = a_. Note that this auxiliary system could be 
either a phase or transmon qubit, as each have a similar 
level structure, in that W12 = wo2 ^ ^oi < "^oi, and A « 
•\/2. In addition, they are both tunable by external flux 
pulses, which we will use in our construction. 

An energy level diagram is shown in Fig. [sjja), where 
we have used the convention to label the system by |g, n) , 
were g = 0, 1, 2 is the state of auxiliary qutrit and n is the 
photon number (or Fock state). This is a generalization 
of the classic Jaynes-Cummings Hamiltonian to a three- 
level artificial atom coupled to a resonator. 

This Hamiltonian above conserves the excitation num- 
ber 



iV = ata+(|i)(i|+2|2)(2|), 



(10) 



so that we can break the problem into an infinite set of 
(up to) 3x3 blocks. Using the notation \q,n) for qubit 
states q — 0,1,2 and Fock states n, the ground state 
|g = 0, n = 0) is unique and set to have zero energy. The 
first-excited subspace of = 0, n = 1) and |g = 1, n = 0) 
is governed by the Hamiltonian 



UJr 9 

9 ^01 



(11) 



The remaining states involve 3x3 matrices for the states 
|0,n), |l,n-l), |2,n-2): 



5V"- 







"^2 = I 9V^ {n — l)oLJr + ujol 9y^\Jn - 1 

9^\Jn — 1 (n — + 10^2 

■(12) 

Assuming that we are away from the avoided crossings 
wqi = w,-, W12 — w,. or W02 — 2ujr, we can apply pertur- 
bation theory to to find the following for the energies 
E ■ 



LiJr — LOqi 



El^n/f^ ~ nUJr + UjQi + {n + 1) - 



jr 

LOqi — UJr 



2\2 



UJr - UJ12 



E2,n/tl ~ nwr + a;o2 + (n -f 1)- 



2\2 



.rA' 



UJ12 — UJr 



(13) 



(14) 
(15) 



The shift of the eigenvalues is the AC Stark shift, and 
have been seen for coupling of a qubit to both quantum 
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FIG. 4: The number-state-dependent Stark shift a;oi,n — <^oi 
as a function of the resonator frequency LUr- The perturba- 
tive Stark shift is shown for n = — >■ 5 as a function of 
the resonator frequency iJr/2-K, for typical qubit parameters 
ojoi/Stt — 7 GIfz, UJ12/2-K = 6.58 GHz, A = 1.46 and coupling 
g/2n = 35 MHz. 
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FIG. 5: (a) Encoding and decoding sequence for resonator 
logic gate. The number-state-dependent qubit rotations R, 



(1) 



("2) 

and Rqi , in the dispersive regime, select out the Fock states 
n = 1 and n = 2, while the R12 transition prepares the state 
for the logic gate, (b) The swap gate S performs the logic 
gate in the resonant regime tJi2 — oJr- 



section), as was recently used for NOON state prepration 



and classical fields. In the dispersive regime, we can ef- 
fectively eliminate state g = 2 to have the modified level 
diagram shown in Fig. [sjb). 

As a consequence of this shift, the transition between 
qubit states depends on the photon number. Defining 



-{2n+l) 



(16) 

This shift of the qubit transition is indicated in Fig.[3j[b). 
This Stark shift can be used to provided a number-state- 
dependent transition, effectively a controlled-rotation of 
the qubit based on the Fock state of the resonator, by 
applying an additional microwave field to the qubit of 
the form "Hdrivc — hfl{a+ + tT_) cos wi. 

The frequency shift Wq"'' — wqi is shown in Fig. jij with 
typical experimental parameters. As described in Koch 
et al. [33], there are three special regions in this figure: 
iOr < LO12, oj\2 < uJr < wqi, and wol < i^r- This is 
different from what would be expected for a resonator 
coupled to a two-level system, which would have only 
two regions, one with positive shift and one with negative 
shift. The middle region with positive shift is known as 
the "straddling" regime, and has the largest value, while 
the negative regions have smaller shifts. The divergences 
in Fig. |4] are at the resonant conditions ujj, = wqi or lo^ = 
0Ji2. These are in fact avoided crossings where the states 
(and transitions between them) are more complicated. 

These avoided crossings can also be used for control. 
By rapidly shifting the qubit frequency to one of these 
anticrossings by a "shift" pulse, swapping between the 
hybridized states occurs ||5j • This has been used to swap 
excitations from the qubit to the oscillator 04] and to 
prepare Fock states and their superpositions [121 HI] ■ We 
will consider the anticrossing at UJ12 = uj^ (see the next 



B. Qudit operation 

Having defined the quantum system, we now illustrate 
how the dispersive and resonant interactions can be used 
for arbitrary single-qudit operation. Consider a quantum 
state: 



IV-o) = |0)® ^c„|n). 



(17) 



n=0 



We begin by performing a rotation between neighbor- 
ing Fock states j and j + 1. To do this, we first ap- 
ply a number-state-dependent 7r-pulse, conditioned on 
the photon state n = j, performing the transformation 

|0,n) — >• \6nj,n); this will be called Rq^ . This is fol- 
lowed by a TT-pulse on the 5 = 1 — > 2 transition, called 
R12, after which another number-state selective tt pulse 
is performed, conditioned on the photon state n = j + 1. 
The net result of these operations is to transform \tpo) 
into 



C^cncodcl'Ao) 

c,|2,j) -f cj+i|l,j + l) + J2 c„|0,n). 



(18) 



where Uo 



ode 



Rli'i^''Ri2Roi- This has selected out 



the |j), \j + 1) subspacc of the resonator; this sequence is 
illustrated in Fig. [sjja) for j = I. 

The system is now configured to the resonant regime, 
with the resonator frequency equal to the (7 = 1 — )• 2 
transition frequency UJ12, as shown in Fig. [sjjb). This 
can be done by dynamically tuning the qubit frequency. 
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and was the key step in the NOON state experiment 
[24] , The subsequent evolution is a two-state oscihation 
between |1, j + 1) and |2, j), so that 

|V2> - ^WIV-i) 

= c,|2,j) +5j+i|l,j + l> + J2 c„|0,n), 

(19) 

where 6 — Xgt and the new amphtudes are 

Cj ~ cos 9cj — ism 9cj^i 
Cj+i = cos0Cj-|_i — isinScj 

To remove the entanglement between the qubit and 
the oscillator, we reverse the encoding step. That is, we 
perform the number-state-dependent 7r-pulse Rqi^^\ the 
q = 2 — 1 transition R12 , and finally R^^^^ . The net result 
is to map \2,j) — > \0,j) and |l,j + l) |0, j + so that 

1^3) = C^dccodclV'2) 

= cj\0,j) +cj+i\0,j + 1) + c„|0,n), 

(21) 

where C/docodc = Rqi RuRqi'^^ ■ This has achieved the 
desired rotation, TZjj^i{9). In short, we have found 

H,^,+m ^ Rg^R,2R^^+'^Si9)Rg-''^R,2R^,\ (22) 

As alluded to above, any Fock-state rotation TZj,k{9) 
can be implemented by using the nearest-neighbor rota- 
tions Tljj+i{9). This is done by swapping state ampli- 
tudes along paths in a "coupling graph" , as described in 
|38j . For example, we can extend our construction to the 
rotations 

7^JJ+2(0) = Rg^Ri2S{7r)Ri2R^^+^^S{9) (23) 
xRq^^^ Ri2S{Tr)Ri2RQi . 

and 

nj,,+3{9) = 4f i?i25(7r)i?i25(7r) (24) 
xR,2Rg'''^Si9)Rg+'^R,2 

X5(^)i?125(^)i?124'l^ 

Note that each of these has the form C/docodc'S'(0)[/oncodc: 
we first transform the state by encoding it into a par- 
ticular set of qudit states (suitably entangled with the 
auxiliary), perform a swap, and then decode the state so 
that the net result is a transformation of the qudit state 
alone. 



IV. NUMERICAL SIMULATION 

We have solved the Schrodinger equation for a four- 
level system coupled to a resonator. The lowest few en- 
ergy levels En, n = 0, 1, 2, • • • , for this system are shown 




FIG. 6: The energy levels Ei and E2 for the coupled 
qubit-resonator system, as a function of the qubit frequency 
ojoi. These levels correspond to the single-excitation sub- 
space, with eigenstates l^*!) « |0, 1) and {'^2) ~ |1,0) for 
UJ01/2-K > 7.3 GHz. 




FIG. 7: The energy levels £"3, £4, and E5 for the coupled 
qubit-resonator system, as a function of the qubit frequency 
LjQi. These levels correspond to the two-excitation subspace, 
with eigenstates 1*3) ^ |0,2), I'I'i) ~ |1,1), and I'I's) ~ |2,0) 
for u}oi/27v > 7.6 GHz. 



in Figs. |6] and [7) We have used four levels for the aux- 
iliary and ten for the resonator, with parameters similar 
to transmon-style qubits: (wqi — Wi2)/27r = 420 MHz, 
(ujQi - W23)/27r = 910 MHz, LOr/2'K = 7 GHz, and 
g/27r — 35 MHz. These are similar to recent experi- 
ments, and the resulting levels are very similar to the 
energy levels for coupled phase qubits US] . 

Three avoided crossings are indicated in Fig. [Tj The 
first has wqi = w^, while the second, W02 = 20;^, is a 
second-order crossing. The gate described above uses 
the third avoided crossing at W12 = Wr- Away from these 
crossings, we can define a Stark shift for transitions be- 
tween states predominantly composed of the uncoupled 
eigenstates \q,n). In the following, we will use these 
"dressed" eigenstates to characterize our logic gate. The 



6 



150r 



100 

X 

1, 50 



0^ 



-50 



-100 





\\\ 

p 











^ CD /2ii (GHz/'^ 



FIG. 8: The number-state-dependent Stark shift ojq"' — ojqi 
as a function of the qubit frequency a;oi. The numerically 
calculated Stark shift is shown for n = — > 7. 



(n) 

Stark shift — wqi as a function of qubit frequency 
for the various Fock states is shown in Fig. [8j the ad- 
ditional structure in the straddling regime is due to the 
second-order crossing. 

To illustrate the logic gate sequence described above, 
we start in the straddling regime with a;oi/27r = 
7.28GHz, and implement the control sequence shown in 
Fig. [21 The Ion ger microwave pulses implement the 
number-state-dependent rotations, the shorter pulses the 
g = 1 — > 2 transition, while the upper shift pulse imple- 
ments the swap operation. All of the microwave pulses 
use a truncated Gaussian profile fTFj . Here the qubit fre- 
quency is shifted from aJoi/27r = 7.28 7.49 GHz and 
back, causing the exchange |1,1) — |0, 2). The ampli- 
tude or the timing of this shift pulse can be adjusted for 
an arbitrary rotation; here we have chosen to perform a 
full swap rt = 0— >-lorn = l— >-0. For this choice of 
system parameters, the complete sequence takes 346ns. 
In terms of the dressed eigenstates, this swap is between 
the energy eigenstates and |\E'2). This sequence can 
be extended to perform swaps between any neighboring 
Fock states with a similar control pulse. 

Solving the time-dependent Schrodinger equation, the 
probabilities pk = |(vt'fc|'I'(t))p (for_the first few eigen- 
states) are shown in Figs. 



10 and 11 



using initial condi- 
tions appropriate to n = and n — I, respectively. The 
swap probability for these states is ~ 0.99 , while the 
higher Fock states are unaffected (with fidelity > 0.95). 
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FIG. 9: Control sequence for resonator logic gate Uo,i- The 
upper line indicates the qubit frequency as a function of time, 
while the lower indicates the microwave pulses required to 
encode and decode the appropriate Fock states. The various 
steps of the sequence are labelled. 
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FIG. 10: Time-dependent probabilities = |(*fc|'I'(t))p for 
the swap n = — ^ 1. The solid black lines are for k = 1 and 
2, while the dashed and gray lines are for = 3 — > 6. In this 
simulation, the initial state 1^(0)) = |>I'i) is evolved using the 
control sequence of Fig. |9] 



state-dependent rotations to encode and decode the Fock 
state I j, k) to be coupled. Starting with 



V. TWO-QUDIT GATE 

An extension of this scheme to two-qudit, and hence 
arbitrary quantum computation, will now be described. 
To perform the two-qudit operation Z/,-,fc, we return to 
the circuit of Fig. [Sj with an auxiliary qutrit for each 
qudit, and the qudits are coupled together. We denote 
the states by \qai qbjna^rib) ■ We again use the number- 



lV'o) = |0,0)®^ \n,m), (25) 

n,m 

the encoding operation [/encode = ^^4 01^5 01 Prepares 
the system in the state 

= C/oncodc|l/'o) = Cj-fc|l,l,J,fc) -I- (26) 
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FIG. 11: Time-dependent probabilities Pk = for 
the swap n = 1 — > 0. The sohd black lines are for k = 1 and 
2, while the dashed and gray lines are for fc = 3 — >■ 6. In this 
simulation, the initial state |^(0)) = |^2) is evolved using the 
control sequence of Fig. |9] 
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FIG. 12: Control sequence for resonator two-qudit logic gate 
Wn,n. The upper two curves indicates the microwave pulses 
and qubit frequency as a function of time for system B, while 
the lower two indicate the qubit frequency for system B. The 
various steps of the sequence are labelled, with C correspond- 
ing the controUed-phase gate operated at wa.oi = ub,i2 (see 
text). 



where 

+ Cn,„,\0,0,n,m). (27) 

This operation has selected out the oscillator states with 
Ua = j and rib = k. 

A quantum logic operation can now be performed on 
the qutrits, specifically a controUed-phase gate C{9) of 
the form \qa,qb) ^ e*'^"''|ga, ^fc), where (f>ab = for 
9a = 9h = 1 and zero otherwise. This generates the 
transformation 

IV2) =C(0)|V'i)=e^'c,-fe|l,l,j-,A;) + |^V>. (28) 
Finally, by using J/dccodo = ^a,oi^b,oi' ^^^^ 

IV'3) = C/dccodo|?A2) 

= e*^Cj,fc|0, 0, j, fc) -I- ^ c„^„|0, 0,n,rn), 

n,m^(j,/c) 

(29) 

returning the encoded states to the resonator. In sum- 
mary, we have shown that 

l^j,k{0) = ^aoi-^b,oiC'(^')^aoi-Rb,oi- (30) 

Combining logic gates of this form with single-qudit op- 
erations allows for universal quantum computation over 
an arbitrary number of qudits [38] 

The controlled-phase gate C{6) between the two 
qutrits can be implemented by shifting their fre- 
quencies so that uiAfii = i^B,i2- The interaction 



figab{'^+<7- + <^-0'+) now leads to the resonant exchange 
1 1,1) -i|0,2) -|1,1) [5]. This TT-phase shift can 
be adjusted to any value by using a nonzero detuning 
'^A.oi ^ ^B,i2 [4S;, or by an adiabatic implementation 
j2j. The full control pulse, assuming gab/'^^'K = 35MHz, is 
shown in Fig. [12] taking a total time of 150 ns. 



VI. DECOHERENCE AND MEASUREMENT 

Resonators have very attractive coherence properties, 
with the potential for more complex qudit operations 
than their qubit counterparts. Resonators have shown 
nearly ideal decoherence dynamics [20 , described mainly 
by energy loss. On-chip resonators typically have co- 
herence times greater than 1 /is, while recent three- 
dimensional cavities have shown qubit coherence times 
Tq greater than 10 /is and resonator coherence times 
greater than 50 /xs [47]. However, the n-th excited state 
of the resonator decays with a rate n/Tj., proportional to 
the Fock state number. A reasonable conclusion is that a 
"good" resonator qudit could have d < 1+Tr/Tq « 6—10. 
We will numerically simulate the gate sequences de- 
scribed above to verify this conclusion. 

Resonator qudits will also require require a means to 
readout the resonator state. The simplest method would 
use the quantum Rabi oscillations for wqi = uJr, as in 
the experiments of Hofheinz et al. [HI |2T]. Here the 
exchange of energy between qubit-resonator states |0,n) 
and |1, n — 1) occurs with (angular) frequency gy/n. This 
allows the populations of the various Fock states to be 
found by collecting a suitably long time-series and Fourier 
analysis. An alternative method would use the number- 
state-dependent Rabi transitions to implement the non- 
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demolition method of Johnson et al. ^3]. A sequence of 
such transitions applied to a qubit initially in its ground 
state would allow the populations of the various Fock 
states to be determined, one by one. Both methods 
would require repeated qubit measurements to estimate 
the Fock state probabilities. 



A. Decoherence Simulation 

We model decoherence using the Lindblad master 
equation 



dp 



1 



(31) 

B 

— 1 



with up to four Lindblad operators L\ ~ a_, L2 ~ u 
L3 = a, L4 = 6, and rates Ai = A2 = l/Tlj and A3 = 
A4 = l/Tr- To simplify the calculation, we transform to 
an interaction picture and keep only the resonant terms 
in H for each step of the logic gate. 

For the single-qudit logic gate, this entails the sequence 
of interactions 



^int/^ 



\^2{\l){2\ 



7t. 



|l>(0|)®|j)01 fori?; 
|2)(1|)®/ 



0) 

01 ' 



for i?i2, (32) 

for s{e), 



where / is the identity operator for the resonator. The 
resulting swap probabilities for the single-qudit rota- 
tions Un,n+i{'^) are shown in Fig. 13 for Fock states 
n = — >■ 7. These simulations use a quantum trajectories 
approach to integrate the master equation, with 1024 tra- 
jectories. The upper curve is for state-of-the-art coher- 
ence times, while the lower is for typical on-chip circuits. 
These results are consistent with a loss of coherence pro- 
portional to e-0-66T/T,g-nT/T,^ where T = 342 - 346ns 

is the total time for the single-qudit gate. As discussed 
above, the resonator Fock states with n < Tr/Tq will have 
errors of about the same order as a single-qubit gate of 
the same duration. 

For the two-qudit controUed-phase gate, the interac- 
tion Hamiltonians are 



'<yx(^l(^\j){j\ 



for R 
for R 



U) 

A, 01' 

(i) 

S,01' 



= \/2.g(|ll)(02| + |02)(ll|)(^/(g)/ for C(6l), 

(33) 

where I a, Ib, and / are the identity operators for the 
auxiliaries A, B, and a resonator, respectively. The 
resulting (worst-case) fidelities for the two-qudit gate 
^n,n(7r) are shown in Fig. W^ again calculated using the 
quantum trajectories method for the master equation. 
These results are somewhat better than the single-qudit 
gate, proportional to e~^/^''e~^"^/'^'', here with a smaller 
overall time T = 160 ns. 




3 4 
Fock State # 



FIG. 13: Swap probabilities for the single-qudit rotations 
i7n,n+i(7r), for n = — > 7. The upper dots are for coherence 
times Tq = 10 ps, Tr — 50 ps, while the lower squares are 
for coherence times Tq — 1 ps, Tr = 10 ps. Other relevant 
parameters are f2i/27r = 6.67MHz, il2/2'K = 25MHz, and 
g/2n — 35MHz. The curves are guides for the eye. 




3 4 
Fock State # 



FIG. 14: Gate fidelities for the two-qudit phase gate (tt) 
for n = — > 7. The upper dots are for coherence times Tq = 
10 ps, Tr — 50 ^s, while the lower squares are for coherence 
times Tq = 1 ps, Tr = 10 ps. Other relevant parameters 
are n/2-K = 6.67MHz and gab/2TT = 35MHz. The curves are 
guides for the eye. 



VII. CONCLUSION 

We have presented an approach to quantum computa- 
tion using the multilevel Hilbert space of a resonator as 
a new type of qudit. This approach is based on resonant 
and dispersive interactions that have been demonstrated 
experimentally, and can be extended to multi-resonator 
logic gates. Thus, a successful demonstration of this 
scheme will open up a number of interesting questions 
in quantum information processing. 

First, while there is a great deal known about the 
theory of quantum circuits and algorithms for qubits 
|27) . much remains to be learned about qudit algorithms. 
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While the asymptotic complexity should be identical [55] , 
these results indicate that an arbitrary unitary gate on 
n qudits requires d" elementary gates (exponential in n) . 
Efficient quantum algorithms, using specific gates such 
as the quantum Fourier transform, can be implemented 
using a polynomial number of qubit gates. An interest- 
ing problem would be to determine if qudit constructions 
for the Fourier transform can be more efficient than the 
qubit constructions. 

Second, our analytical approach leaves open questions 
about optimization of operations in this larger Hilbcrt 
space. There may be interesting approaches to construct 
a given unitary of interest that is more efficient than the 
two-level reduction used here. In particular, the number- 
state-dependent transitions dominate the operation time. 
While this time can likely be decreased by using larger 
couplings or more sophisticated microwave pulses, per- 
haps using multiple frequencies |15^ or multiple quadra- 
tures [16 , other approaches may be necessary. Optimal 
control methods [IS] for this system may lead to such 



alternative approaches, and would be important for op- 
erations in the presence of decoherence. 

Finally, the two measurement approaches presented 
both use qubit measurements to read out the resonator 
states. Either approach should allow for full tomogra- 
phy of the resonator logic gates, at the expense of having 
to perform a large number of experiments to determine 
the state of the resonator. An open question is how to 
extract this information in the most direct and efficient 
manner. These and other issues will be fruitful tests of 
our understanding of quantum control and measurement 
of multilevel quantum systems. 
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